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ABSTRACT 

We present the analysis of 12 high-resolution galactic rotation curves from The HI Nearby Galaxy Survey (THINGS) in the context 
of modified Newtonian dynamics (MOND). These rotation curves were selected to be the most reliable for mass modelling, and they 
are the highest quality rotation curves currently available for a sample of galaxies spanning a wide range of luminosities. We fit the 
rotation curves with the "simple" and "standard" interpolating functions of MOND, and we find that the "simple" function yields 
better results. We also redetermine the value of «o, and find a median value very close to the one determined in previous studies, ao 
= (1.22 ± 0.33) x 10~ 8 cm s~ 2 . Leaving the distance as a free parameter within the uncertainty of its best independently determined 
value leads to excellent quality fits for 75% of the sample. Among the three exceptions, two are also known to give relatively poor fits 
also in Newtonian dynamics plus dark matter. The remaining case (NGC 3198), presents some tension between the observations and 
the MOND fit, which might however be explained by the presence of non-circular motions, by a small distance, or by a value of ao 
at the lower end of our best-fit interval, 0.9 x 10~ 8 cm s~ 2 . The best-fit stellar M/L ratios are generally in remarkable agreement with 
the predictions of stellar population synthesis models. We also show that the narrow range of gravitational accelerations found to be 
generated by dark matter in galaxies is consistent with the narrow range of additional gravity predicted by MOND. 
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1. Introduction 

The current dominant paradigm is that galaxies are embedded 
in halos of cold dark matter (CDM), made of non-baryonic 
weakly-interacting massive particles (e.g., Bertone et al. 2005). 
However, an alternative way to explain the observed rotation 
curves of galaxies is the postulate of Milgrom (1983) that for 
gravitational accelerations below a certain value ao, the true 
gravitational attraction g approaches (gN^o) 1 ^ 2 where g^ is the 
usual Newtonian gravitational field (as calculated from the ob- 
served distribution of visible matter): this paradigm is known as 
modified Newtonian dynamics (MOND). 

MOND explains successfully many phenomena in galax- 
ies, among which the following non-exhaustive list: (i) it pre- 
dicted the shape of rotation curves of low surface-brightness 
(LSB) galaxies before any of them had ever been measured (e.g. 
McGaugh & de Blok 1998); (ii) tidal dwarf galaxies (TDG), 
which should be devoid of collisionless dark matter, still exhibit 
a mass-discrepancy in Newtonian dynamics, which is perfectly 
explained by MOND (Gentile et al. 2007); (iii) the baryonic 
Tully-Fisher relation (e.g., McGaugh 2005), one of the tightest 
observed relations in astrophysics, is a natural consequence of 
MOND, both for its slope and its zero-point; (iv) the first realis- 
tic simulations of galaxy merging in MOND were recently car- 
ried out, notably reproducing the morphology of the Antennae 
galaxies (Tiret & Combes 2008); (v) it naturally explains the 
universality of "dark" and baryonic surface densities within one 
core radius in galaxies (Donato et al. 2009, Gentile et al. 2009). 

Recent theoretical developments have also added plau- 
sibility to the case for MOND through the work of, e.g., 
Bekenstein (2004), Sanders (2005), Zlosnik, Ferreira & 



Starkman (2007), Halle, Zhao & Li (2008), and Blanchet & 
Le Tiec (2008), who have all presented Lorentz-covariant the- 
ories yielding a MOND behavior in the weak field limit. 
Although still fine-tuned and far from a fundamental theory 
explaining the MOND paradigm, these effective theories re- 
markably allow for new predictions regarding cosmology (e.g., 
Skordis et al. 2006) and gravitational lensing (e.g., Angus et 
al. 2007, Shan et al. 2008). For reviews of MOND's suc- 
cesses and weaknesses, both at the observational and the- 
oretical level, as well as comparisons with dark matter re- 
sults, see McGaugh & de Blok (1998), de Blok & McGaugh 
(1998), Sanders & McGaugh (2002), Bruneton & Esposito- 
Farese (2008), Milgrom (2008), Skordis (2009), Famaey & 
Bruneton (2009), Ferreira & Starkman (2009). 

One thing the MOND paradigm does not directly predict, 
though, is the shape of the interpolation between the MONDian 
regime where g <K ao and the Newtonian regime where g » ao, 
as well as the actual value of the acceleration constant ao- The 
latter is in principle a free parameter, but once its value has been 
determined by some means, it must be identical for every astro- 
nomical object. Large variations of ao would invalidate MOND 
as a fundamental paradigm underpinned by new physics. Let us 
note that, as shown in Begeman, Broeils & Sanders (1991) fits 
with variable ao and fixed distance D are essentially identical 
to fits with fixed ao and variable D because the observed total 
gravitational acceleration is proportional to l/D. Ideally, the fit- 
ted distance should however generally conform to the indepen- 
dently determined one (e.g., Cepheids-based or RGB tip-based). 
Finally, a consequence of the absence of galactic dark matter 
within the MOND context is that the dynamical mass-to-light 
ratio that is derived from a rotation curve fit should agree with 
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the true stellar mass-to-light ratio of the stellar disk (and some- 
times bulge), as inferred from e.g. observed colours and stellar 
population synthesis models. 

Here, we use results from The HI Nearby Galaxy Survey 
(THINGS; Walter et al 2008), which consists of high-resolution 
HI observations of a sample of 34 nearby galaxies, in order 
to constrain the transition function of MOND. In particular, 
we show that some individual galaxies that had been claimed 
to be potentially problematic for MOND such as NGC 2841 
(Begeman et al. 1991) can yield good fits with the "simple" in- 
terpolating function. 

We use a subset of the THINGS galaxies for which rotation 
curves could be derived in de Blok et al. (2008), restricting our- 
selves to galaxies which are not (obviously) dominated by non- 
circular motions. In Sect. 2, we summarize the popular choices 
that have been proposed in the literature for the transition be- 
tween the MONDian and Newtonian regimes, in order to con- 
front these different transitions with THINGS rotation curves. In 
Sect. 3, we explain how we selected the subsample of galaxies 
that we model in the context of MOND. Sect. 4.1. then presents 
the results for the value of the acceleration constant ao, while 
Sect. 4.2. and 4.3. present the comparison of the rotation curve 
fits for the different transitions, especially the best-fit mass-to- 
light ratios and distances. Finally, in Sect. 4.4., we discuss NGC 
3198, the only cases where the MOND fits perform significantly 
worse than dark matter fits in the context of Newtonian dynam- 
ics. Conclusions are drawn in Sect. 5. 



2. MOND and its interpolating function 

The MOND paradigm stipulates that the Newtonian acceleration 
gN produced by the visible matter is linked to the true gravita- 
tional acceleration g by means of an interpolating function //: 

4-)g = SN, (1) 

where //(x) ~ x for x <SC 1 and //(x) ~ 1 for x » 1 (and g = \g\). 
However, this expression cannot be exact for all orbits and all ge- 
ometries, since it does not respect usual conservation laws. Such 
a modification of Newtonian dynamics could come at the classi- 
cal (non-covariant) level from a modification of either the kinetic 
part or the gravitational part of the Newtonian action (with usual 
notations; 0^ being the Newtonian gravitational potential): 

S = J ^pv 2 d 3 xdt - J {p(f> N + ^^^jd 3 xdt, (2) 

where modifying the first term is referred to as modified in- 
ertia and modifying the second term as modified gravity. 
Milgrom (1994) has shown that within the modified inertia 
framework, Eq. 1 was exact only for circular orbits (for other 
orbits, predictions are difficult to make since the theory is non- 
local). On the other hand, Bekenstein & Milgrom (1984) have 
shown that within a modified gravity framework where |V0Af I 2 
is replaced by a?F(|V0| 2 /a?) in Eq. 2 (<p being the MONDian 
potential and F' = //), the right-hand side of Eq. 1 had to be 
replaced by g^ + s where s is a solenoidal vector field deter- 
mined by the condition that g can be expressed as the gradi- 
ent of a MONDian potential. Milgrom (2010) has proposed an- 
other modified gravity formulation in which |V0^| 2 is replaced 
by 2V0 ■ V(f> N - alQ(\V(f> N \ 2 /al) (tp being the MOND potential, 
<Pn remaining the Newtonian one, and l/Q' = //): in this case, 
the solenoidal field to be added to the right-hand-side of Eq. 1 is 



different from the one in the Bekenstein & Milgrom formulation 
(see also Zhao & Famaey 2010). 

Although Brada & Milgrom (1995), Famaey et al. (2007) and 
Zhao & Famaey (2010) have shown that the expected differences 
in the predictions of the various formulations for rotation curves 
are not very large, they can be of the same order of magnitude as 
the differences produced by different choices for the //-function. 
In order to constrain // within the modified gravity framework, 
one should calculate predictions of the modified Poisson formu- 
lations of Bekenstein & Milgrom (1984) or Milgrom (2010) nu- 
merically for each galaxy model, and for each choice of parame- 
ters. This is left for further works and we choose here to concen- 
trate on the modified inertia formulation for circular orbits given 
by Eq. 1 . 

It is worth noting that other interesting constraints on MOND 
and its //-function could come from studies of the effect of the 
galactic gravitational field on the dynamics of the inner Solar 
System (Milgrom 2009), or from studies of the dynamics per- 
pendicular to the galactic disk at the solar position (Bienayme et 
al. 2009). 

Various choices for the shape of the //-function have 
been proposed in the literature (see especially Milgrom & 
Sanders 2008 and McGaugh 2008), but we rather concentrate 
here on the two most popular choices that have been studied 
so far. The "standard" //-function (Milgrom 1983) yields a rela- 
tively sharp transition from the MONDian (x <k 1, where x - 
g/ao and g is the gravitational acceleration) to the Newtonian 
(x » 1) regime: 

Kx) = , X , . (3) 

Vl + X 

while the "simple" //-function (Famaey & Binney 2005; Zhao & 
Famaey 2006) yields a more gradual transition: 

1 + x 

Fig.[T]displays those two //-functions as a function of x. Let 
us note that the simple function predicts that a constant acceler- 
ation equal to «o nas to be added to the Newtonian gravitational 
acceleration for g » ao. This is, for the values of ao compatible 
with galaxy rotation curves (Sect. 4.1), in strong disagreement 
with orbits of planets in the inner Solar System, and especially 
with measures of the perihelion precession of Mercury. A solu- 
tion is to use an "improved simple" //-function that rapidly inter- 
polates between the simple and standard ones for values of the 
gravitational acceleration g > 10oo (i-e. a higher value than those 
g which are probed by galaxy rotation curves). Such an improved 
simple function is shown as an example on Fig. |T| Nevertheless, 
we use the standard and simple //-functions hereafter (keeping 
in mind that the latter should be modified in the strong gravity 
regime) of Eqs. 3 and 4 in order to perform our MOND fits to 
THINGS galaxy rotation curves. 

3. The sample 

We use a subset of galaxies in the THINGS survey for which 
rotation curves were derived in de Blok et al. (2008). We restrict 
ourselves to galaxies which are not (obviously) dominated by 
non-circular motions. This means we omit the bright disk galax- 
ies NGC 3031 and NGC 4736. While not necessarily dominated 
by non-circular motions, we also omit NGC 2366, IC 2574 and 
NGC 925. These galaxies have a neutral gas distribution that is 
dominated by holes and shells, the signature of which remains 
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Fig. 1. The interpolating //-functions of Eq. 3 (standard, dashed 
line) and Eq. 4 (simple, dotted line) are displayed as a function 
of x = g/ao. An improved simple function (solid line) interpolat- 
ing between simple and standard for x > 10 is also presented in 
order to show that a transition behaviour governed by the simple 
//-function in galaxies (where x < 10) can a priori be in accor- 
dance with the Solar System constraints (where x » 10). 
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Fig. 2. Best-fit ag values (using the simple interpolating func- 
tion) vs. central surface brightness in the 3.6 /mi band. 



visible in the radial profile of the neutral gas. In these dwarf 
galaxies, the neutral gas profile dominates the total radial bary- 
onic mass distribution, and as the MOND prediction is derived 
from the observed radial surface density distribution these re- 
maining signatures of the holes and shells could possibly lead to 
erroneous results. The analysis of these specific rotation curves, 
interesting as they might be, is left for a forthcoming paper. 

This leaves us with a total of 12 galaxies: some of these have 
already been discussed in Bottema et al. (2002), but with the 
higher resolution data available, and constrained stellar mass- 
to-light ratios as observed in the Spitzer IRAC 3.6//m band, we 
are able to perform a slightly more stringent test for MOND. 
In fact, the 12 rotation curves that we use here are the high- 
est quality rotation curves currently available (in terms of spa- 
tial/spectral resolution and extent) for a sizeable sample of galax- 
ies spanning a wide range of luminosities, and they therefore 



Table 1. Galaxy distances and the methods used to determine 
them. References are: 1: Vinko et al. (2006). 2: Macri et al. 
(2001). 3: Drozdovsky & Karachentsev (2000). 4: Karachentsev 
et al. (2002). 5: Kelson et al. (1999), Freedman et al. (2001). 6: 
Walter et al. (2008). 7: Rawson et al. (1997). 8: Karachentsev et 
al. (2004). 9: Hughes et al. (1998). 



Name 


distance(Mpc) 


Method 


Ref 


NGC 2403 


3.47 = 


t0.29 


SN, Cepheids 


1 


NGC 2841 


14.1 = 


t 1.5 


Cepheids 


2 


NGC 2903 


8.9 + 


2.2 


brightest stars 


3 


NGC 2976 


3.56 = 


t0.36 


tip of the RGB 


4 


NGC 3198 


13.8 = 


t 1.5 


Cepheids 


5 


NGC 3521 


10.7 = 


= 3.2 


Hubble flow 


6 


NGC 3621 


6.64 = 


= 0.70 


Cepheids 


7 


DDO 154 


4.30 = 


t 1.07 


brightest stars 


8 


NGC 5055 


10.1 = 


1=3.0 


Hubble flow 


6 


NGC 6946 


5.9 ± 


1.5 


brightest stars 


8 


NGC 7331 


14.72 


± 1.29 


Cepheids 


9 


NGC 7793 


3.91 = 


= 0.39 


tip of the RGB 


8 



Notes. NGC 2403 is the only galaxy whose distance differs from Walter 
et al. (2008). We used the estimate of the distance given in Vinko et al. 
(2006) because it comes from more numerous and more recent data. 



represent an important test for MOND or any theory that aims 
at fitting galaxy kinematics. In de Blok et al. (2008) some differ- 
ences between their rotation curves and those of previous publi- 
cations were highlighted. These differences could be caused by 
the different approach taken by de Blok et al. (2008) to derive 
the velocity field: they fit the velocity profiles using third-order 
Gauss-Hermite polynomials, instead of the more conventional 
intensity-weighted mean. 

The two baryonic contributions to the rotation curve, nec- 
essary to compute the MOND mass model, were derived by de 
Blok et al. (2008) as follows. First, the shape (but not the ampli- 
tude) of V s t a r S , the contribution of the stars to the rotation curve, 
was derived from the observed 3.6 //m surface brightness profile, 
and slightly modified to account for the observed (J - K) colour 
gradients as a function of radius (which are an indication of a 
radially varying stellar mass-to-light M/L ratio). Although there 
might be some contamination due to young stars and hot dust, 
this contamination is thought to be a negligible contribution to 
the flux at 3.6/mi, see e.g. Pahre et al. 2004, Li et al. 2007, hence 
the 3.6/mi emission is considered as good a tracer of stellar mass 
as the more commonly used K-band (Zhu et al. 2010). For the 
vertical distribution of the stellar disk, de Blok et al. (2008) as- 
sumed a sech 2 distribution with a scale height of zo = h/5, where 
h is the radial exponential scale length. The amplitude of V s tars 
is scaled according to the global stellar mass-to-light M/L ratio, 
which is left as a free parameter and then compared to the pre- 
dictions of stellar population synthesis models (e.g. Bell & de 
Jong 2001). The contribution of the gaseous disk to the rotation 
curve, V gas , was derived from the observed HI surface density 
profiles, and then corrected for primordial He. 

The galaxy distances are determined by various methods 
(Cepheids, tip of the Red Giant Branch, Hubble flow, bright- 
est stars), and their quoted uncertainties are often close to 10%. 
They are listed in Table [T] 
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4. Results 

4. 1 . The acceleration constant a 

The acceleration constant ao of MOND, though unknown from 
first principles, must be the same for all galaxies, therefore it 
has to be determined empirically, e.g. by fitting rotation curves. 
Begeman, Broeils & Sanders (1991) determined the value of the 
acceleration constant ao to be 1.21 x 1CT 8 cm s~ 2 from mass 
modelling of a number of nearby galaxies with the standard 
//-function of Eq. 3. This value was confirmed by Sanders & 
Verheijen (1998) using a sample of rotation curves of galaxies 
belonging to the Ursa Major galaxy group. However, Bottema et 
al. (2002) noted that using an updated value of the distance to 
the Ursa Major group would bring the value of oq down to 0.9 x 
10~ 8 cm s~ 2 . 

The first fits that we performed were those with ao as a free 
parameter (the stellar M/L ratio being the other free parame- 
ter). The distance in these fits was fixed at the values given in 
Table 1, the most accurate for each galaxy to date, to the best 
of our knowledge. We remind the reader that a fit with ao free 
and the distance fixed is equivalent to a fit with the distance free 
and ao fixed (Begeman et al. 1991), because the observed total 
gravitational acceleration is proportional to 1 /D, where D is the 
distance. 

Swaters, Sanders & McGaugh (2010) find a weak correla- 
tion between the R-band central surface brightness and the best- 
fit value of ao as a result of making MOND mass models of 27 
dwarf and low surface brightness galaxies. They find that lower 
surface brightness galaxies have a tendency to have lower ao. 
In Fig. [2] we look for a similar relation using our best-fit val- 
ues of ao from the MOND fits using the simple //-function. We 
do not find the same correlation: indeed, the best-fit values of 
ao are scattered around the median value without any obvious 
correlation with central surface brightness. A thorough interpre- 
tation of the correlation (or lack thereof) between best-fit ao and 
central surface brightness goes beyond the aim of this paper, but 
our finding might not invalidate the interpretation by Swaters et 
al. (2010) that low surface brightness galaxies could be biased 
towards lower values of ao because of the external field effect 
(e.g. Milgrom 1983). In our sample, apart from DDO 154, the 
galaxies have a relatively high surface brightness. 

The result is that the median values are (1.27 ± 0.30) x 10" 8 
cm s~ 2 for the standard // function (eq. [5]) and (1.22 ± 0.33) x 
10~ 8 cm s~ 2 for the simple // function (eq. ^J, values that are re- 
markably similar to the estimates made in previous studies (the 
uncertainties are calculated following Miiller 2000). However, 
for consistency we will now use these new values in the remain- 
der of the paper. We note that our estimates of ao lie between 
cHo/(2n) ~ 1.1 x 10 -8 cm s~ 2 (where H is the Hubble constant) 
and c VA/(27r) « 1.5 x 10~ 8 cm s~ 2 (where A is the cosmolog- 



ical constant). However, the estimate of oq given by Bottema et 
al. (2002), 0.9 x 10~ 8 cm s~ 2 , cannot be excluded by the present 
data, see Section|4.3| 



4.2. Mass-to-light ratios 

Starting from fixed distances to the 12 galaxies, we found in 
Section 4.1a common median value of a ( > corresponding to each 



interpolating //-function. Using this we perform 6 different types 
of fits to each galaxy rotation curve. For each of the two //- 
functions, we make fits with a fixed value of a and (i) a fixed 
distance, (ii) a distance constrained to lie within the error bars 
from its independent determination, and (iii) a free distance with 




Fig. 3. Stellar M/L ratio in the 3. 6//m band vs. ( J-K) colour. The 
full circles are the results of the MOND fits (using the simple //- 
function of Eq. 4 and ao = 1-22 x 10" 8 cm s" 2 ) with the distance 
constrained within the uncertainties of its independently deter- 
mined value, whereas the solid line represents the predictions of 
stellar population synthesis models (see text for details). 




0.4 0.6 
(B-V) 



Fig. 4. Stellar M/L ratio in the B band vs. (B - V) colour. See 
Fig. [3] for the explanation of line and symbols. 



no constraints. In all cases, the stellar mass-to-light ratio of the 
disk (and bulge if present) is left as a free parameter. All the 
results are listed in Table 2. 

From the x 1 values Q the lack of systematic deviations and 
the small number of highly discrepant data points, one can con- 
clude that the fits are generally good, with only a few exceptions 
(cf. next Section). However, one has to check that, when a bulge 
is present, the mass-to-light ratio of the disk is smaller than that 
of the bulge, and that the stellar M/L ratios are realistic. 

There are five galaxies with a bulge. Using the standard ju 
function, two galaxies have the best-fit stellar M/L of the disk 
larger than the one of the bulge. Of the remaining three cases, 
one is undetermined (NGC 2903: its best-fit M/L of the bulge 
is zero but values larger than the best-fit M/L of the disk give 
almost equally good fits; note also the likely presence of a bar, 
Leroy et al. 2009), and two are realistic (the M/L of the disk 
is smaller than M/L of the bulge). On the other hand, when us- 
ing the simple // function no such problems arise, and the M/L 
of the bulge is always realistic. We therefore conclude, in line 
with Sanders & Noordermeer (2007), that the simple // function 
gives superior fits. For the rest of the paper we will thus only 



1 We note that in de Blok et al. (2008) the velocity difference between 
the approaching and receding side was considered in the error budget of 
each point of the rotation curve. In many cases this results in errorbars 
that are larger than the point-to-point scatter, which implies that the \ 2 
values cannot be used as probability indicators. 
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use the simple // function. This justifies our proposal in Fig. [TJ 
of a fi function that resembles the simple one at typical galac- 
tic gravitational accelerations (and the standard one for higher 
accelerations representative of, e.g, the Solar System). 

In order to check how realistic the fitted stellar M/L ratios 
are, we compared them with the results of stellar population syn- 
thesis models. In Fig. [3] we plot the best-fit "global" M/L ratio 
in the 3.6//m band vs. (J - K): 



M/L = 



(M/ L)diskAiisk + (M/ i)bulgeibulge 
^disk + ^bulge 



(5) 



The solid line represents the population synthesis models predic- 
tion (from Bell & de Jong 2001, using also eq. 4 of de Blok et 
al. 2008) with a "diet-Salpeter" IMF. Although with some scat- 
ter, the points lie close to the prediction, and the M/L ratios in 
the 3.6/mi band vary very gradually with colour, staying con- 
stant around 0.5-1. It is also interesting to compare our results 
in a band where the predicted stellar M/L varies more rapidly 
with colour. To achieve this, we converted our M/L ratios to B- 
band (using the B-band luminosity given in de Blok et al. 2008), 
and we made use of the corrected (B - V) colours given in the 
HyperLeda database (Paturel et al. 2003). For NGC 7793, the 
corrected (B - V) colour was not available and we used the ef- 
fective one. The results are shown in Fig. |4] The best-fit stellar 
M/L of MOND closely follows the predictions of Bell & de Jong 
(2001), in that redder galaxies are best fitted with a higher stellar 
M/L ratio. 

The only three galaxies where the best-fit disk M/L ratio dif- 
fers from the population synthesis one by more than a factor of 
two are NGC 2903, NGC 2976, and NGC 7331 (see Fig. [3] and 
Table 2). In NGC 2903 the MOND fit significantly overpredicts 
the M/L, a phenomenon that was also observed in Newtonian 
mass models with dark matter (de Blok et al. 2008). Let us note 
that, if one would follow population synthesis predictions, the 
predicted stellar disk would be surprisingly very sub-maximum 
for a massive galaxy with a rapidly increasing rotation curve 
with maximum velocity ~ 215 km s . In addition, in the cen- 
tral parts of NGC 2903 there is evidence for a bar (e.g. Leroy 
et al. 2009): the non-circular motions associated to it further 
increase the uncertainties on the mass modelling results (see 
also Sellwood & Zanmar Sanchez 2010). For NGC 2976, on the 
other hand, the MOND fits underpredict the stellar M/L ratio; 
but again, this is also observed in the Newtonian mass models 
with dark matter. In addition, having a stellar disk that strongly 
dominates the kinematics over most of the extent of the gaseous 
disk in a galaxy with maximum velocity ~ 85 km s _1 would also 
be surprising. Mass model degeneracies in this galaxy (in par- 
ticular the MOND mass models with distance free and distance 
constrained) are complicated by the very similar shapes of V stars 
and V gas . In NGC 2976 too, Leroy et al. (2009) find an indication 
for a weak bar. Also in NGC 7331 the MOND fits give a lower 
M/L compared to the expectations from the colours; this is the 
case also in the dark matter fits (see de Blok et al. 2008 where 
a strong dust ring is suggested as a possible explanation for the 
inflated stellar M/L ratios predicted from the colours). 

4.3. Distances 

One then also has to check that the fits with the distances con- 
strained to lie within the error bars from their independent de- 
terminations are of good quality. When this is not the case, it 
means that MOND would predict another distance than what has 
been measured to date. NGC 2841 is, e.g., the most well-known 



and most persistently problematic galaxy for MOND. Begeman 
et al. (1991) pointed out that a good MOND fit could only be 
obtained if the galaxy was a factor ~ 2 further away than the 
Hubble distance of ~ 9.5 Mpc. This large discrepancy was alle- 
viated somewhat when HST Cepheids measurements suggested 
a distance of 14.1 Mpc (Macri et al. 2001), but the discrepancy 
remained. However, these fits were performed with the standard 
//-function (Eq. 3) and not the simple one (Eq. 4). Our fits here 
show that the problem of the distance is solved when using the 
simple //-function (see the reduced \ 2 in Table 2), and that the 
stellar mass-to-light ratio is also in accordance with population 
synthesis models. 

All the fits of the 12 high-quality rotation curves, using the 
simple //-function of Eq. 4, = 1.22 x 10~ 8 cm s~ 2 , and a 
distance lying within the error bars coming from an indepen- 
dent distance determination, are shown in Fig. [5] The fits are 
clearly very good for 9 galaxies (including NGC 3521, whose 
high reduced x 1 value in Table 2 is dominated by the innermost 
two points, which have highly uncertain position angle and in- 
clination, see de Blok et al. 2008). Among the galaxies with 
the 3 least good fits (NGC 3198, NGC 7793, and NGC 2976), 
we do not discuss further NGC 7793 and NGC 2976, since the 
MOND fits present the same failures as the dark matter fits (in 
Newtonian dynamics), therefore we do not consider them as ev- 
idence against MOND. We just briefly note that in NGC 7793 
the value of the inclination angle fitted by de Blok et al. (2008) 
is low and presents large variations in adjacent radii, which re- 
sults in a poorly constrained rotation curve; in NGC 2976 the 
amplitude of the non-circular motions (Trachternach et al. 2008) 
is correlated with the amplitude of the fit residuals. 

Before proceeding with a detailed analysis of possible prob- 
lems with the rotation curve of NGC 3 198 in Sect. 4.4 we finally 
consider the possibility that the true value of aq for all galaxies 



is actually at the lower end of our best-fit interval of Sect. 4. 1 



i.e. a value compatible with the one determined by Bottema et 
al. (2002). As a matter of fact, a good reason for this is that 
the Ursa Major (UMa) galaxy group (e.g., Sanders & Verheijen 
1998, Gentile, Zhao & Famaey 2008) is nowadays thought to be 
at a distance of 18.6 Mpc (Tully & Pierce 2000), implying a best- 
fit value of flo close to the one of Bottema et al. (2002), see e.g. 
Gentile et al. (2008). To get as good fits to the rotation curves 
of UMa galaxies as those obtained with «o = 0.9 x 10~ 8 cm s~ 2 
with a higher value of the order of Oq = 1.2 X 10~ 8 cm s~ 2 , the 
distance of the group should be smaller, of the order of 15 Mpc 
(as originally assumed by Sanders & Verheijen 1998). For this 
reason, we plot in Fig. |6]the fits of the 12 rotation curves using 
«o = 0.9 x 10~ 8 cm s ^and still the simple //-function and the 
distance constrained to lie within the error bars of Table fTJ. As 
can be seen, the fits remain of approximately the same qualit 
apart for 3 galaxies: NGC 2841 and NGC 2403 have worse fits 



4.4. NGC 3198 

An excellent fit can be found with «o = 1-2 x 10~ 8 



cm s 2 by 



leaving the distance unconstrained. As already noted in previ- 



2 In this case the less well fitted galaxy would be NGC 2403. Possible 
effects could be the fact that the r, J, H and K photometric profiles have 
a different shape from the 3.6//m band (see Kent 1987, Fraternali et al. 
2002, de Blok et al. 2008). In addition, Fig. 7 of de Blok et al. (2008) 
shows that the outer parts of the rotation curve are quite uncertain. In 
that paper, the stellar component is also modelled with two separate 
components: the use of two different disks for the stellar contribution 
does not change significanly the results., but the quality of the fit of 
NGC 3198 improves. 
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Fig. 5. Rotation curve fits with the distance "constrained" and «o = 1-22 x 10~ 8 (fi simple d constr. in Table 2). Dashed, dotted, 
and long-dashed lines represent the Newtonian contributions of the gaseous disk, stellar disk, and bulge, respectively. The MOND 
best-fit model is shown as a solid red line. 



ous studies (e.g. Bottema et al. 2002), MOND prefers a smaller 
distance. Fig. [7] shows a very good fit with a distance of 8.6 
Mpc, which is significantly lower than the Cepheids-based one 
of 13.8 Mpc. The distance that one would get by fixing the stel- 
lar M/L ratio to the population synthesis value is also lower 
than the Cepheids value (9.6 Mpc). We note that other meth- 
ods (the Tully-Fisher distance or the Hubble flow distance, both 
less accurate than the Cepheids) also yield lower values (see van 
Albada et al. 1985 and Bottema et al. 2002) than the Cepheids. 
It is however also interesting to note that the regions where the 
fit with the distance constrained is most discrepant (roughly, the 
first and last thirds of the rotation curve) are also those where 
the amplitude of the non-circular motions is higher, taken from 
Trachternach et al. (2008) and they are about of the same order 



of magnitude. The increase of the non-circular motions from the 
middle part of the rotation curve to the outer parts at ~ 10 km s _1 
is also noted in Sellwood & Zanmar Sanchez (2010). The inter- 
pretation is further complicated by the fact that the IRAC 3.6 /im 
image shows what seems to be an end-on bar in the very inner 
part of NGC 3198. The main spiral arms of NGC 3198 originate 
here, and it is not clear how these affect the dynamics. The use 
of two different disks for the stellar contribution (see de Blok 
et al. 2008) does not change significanly the results. We thus 
conclude that a full MOND modelling of the two-dimensional 
velocity field of this galaxy, taking into account the full mod- 
ified Poisson equation(s) of Bekenstein & Milgrom (1984) or 
Milgrom (2010), as well as the non-axisymmetry of the galac- 
tic potential, would lead to a benchmark test for the viability of 
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Fig. 6. Rotation curve fits with the distance constrained to lie within the error bars of Table[T] and ciq = 0.9 x 10 8 cm s 2 . The lines 
are described in Fig. [5] 



MOND as a modification of gravity. We note, however, that a 
lower value of ao would improve the fit of NGC 3198 although 
the fits of NGC 2841 and NGC 2403 would get worse. With 
ao = 0.9 x 10~ 8 cm s~ 2 , a fit equivalent to the one of Fig. [7] 
would be obtained for a distance of 10.3 Mpc instead of 8.6 Mpc, 
closer to the Cepheid-based distance of 13.8 Mpc. Actually, it 
has been mentioned in the literature that there might be a redden- 
ing problem in determining the distance. In Macri et al. (2001), 
the H-band Cepheids distance of NGC 3198 is 11.2 Mpc. So, 
with ao = 0.9 x 10" 8 cm s~ 2 , there might be no problem with the 
rotation curve of NGC 3 198 (see Fig. [8j. 



4.5. Scaling relations 

An interesting way to interpret the above MOND results is to 
phrase them in terms of the usual dark matter framework, con- 
sidering MOND as a phenomenological, empirical, law encom- 
passing the behavior of dark matter in galaxies. The additional 
gravity generated by MOND, compared to the Newtonian case, 
can indeed be attributed to what one would call dark matter in the 
Newtonian context, and this effective matter is called "phantom 
dark matter". For the gravity generated by baryons, we hereafter 
use the mass-to-light ratios from the fits made with the simple 
interpolation function and the distance constrained. 

In Fig.|9| we display the scaling relation known as the Mass 
Discrepancy-Acceleration relation (McGaugh 2004), showing 
that the ratio of enclosed total dynamical mass (in Newtonian 
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Fig. 7. Rotation curve fit of NGC 3198 with a = 1.2 x 1(T 8 cm 



and samples the yu-function. Since for the considered range of 
such gravitational accelerations, the range of "phantom dark 
matter" gravities is not very large, it gives the illusion of a dark 
matter gravity which is more or less constant with radius. 



5. Conclusion 

We re-analysed the ability of the Modified Newtonian Dynamics 
(MOND; Milgrom 1983) paradigm to fit galaxy rotation curves, 
using the most up-to-date high-resolution HI data for nearby 
(d < 15 Mpc) galaxies from the recent THINGS survey (Walter 
et al. 2008). We selected a subset of 12 galaxies not obviously 
dominated by non-circular motions, and yielding the most reli- 
able mass models. 

First, we redetermined the value of the acceleration parame- 
ter in MOND (ao), which is unknown a priori but has to be the 
same for all galaxies. This was done for both commonly used 
interpolating functions /j of MOND. We find a median value of 
a = (1.27 ± 0.30) X 10 " 



8 cm s 2 for the "standard" fi function 



s 1 the distance as a free parameter (p simple d free in Table 2). (Eq. 3), and a = (1.22 + 0.33) x 10~ 8 cm s -2 for the "simple" fi 



The distance is 8.6 Mpc and the stellar mass-to-light ratio in the 
3.6/im band is 1.01. The lines are described in Fig. [5] 



NGC 3198 
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Fig. 8. Rotation curve fit of NGC 3198 withon = 0.9 x 10" 8 cm 



s , a distance of 11.2 Mpc (see Section |4.4[ ), a best-fit stellar 
M/L ratio in the 3.6//m band of 0.76, and the simple fi function. 
The lines are described in Fig. [5] 



gravity) w.r.t. enclosed baryonic mass at any radius is a function 
of the gravity generated by the baryons at this radius. This re- 
lation precisely traces the ju-function of MOND, and the small 
scatter around the line indicate the deviation from the MOND 
behavior, mostly consistent with observational errors. 

Then, following Walker et al. (2010), we plot the gravity of 
(phantom) dark matter as a function of radius for the 12 galaxies 



of the sample. This is plotted on Fig. 10 We remarkably find that 



the additional gravity predicted by MOND is in accordance with 
the mean and scatter in Fig. 1 of Walker et al. (2010). However, 
if one plots (Fig. 11 1 this additional gravity as a function of the 



function (Eq. 4), very close to the value that had been determined 
in previous studies (e.g. Begeman et al. 1991). 

Then, fixing these values for ao, we performed three fits for 
each jj function: with the distance fixed at the value determined 
in an independent way, then by leaving the distance free but con- 
strained within the uncertainties of this distance determination, 
and then with the distance as a free parameter with no constraints 
(Table 2). We find that the MOND fits with the distance "con- 
strained" are of very good quality (Fig. [5}, with three excep- 
tions: two of these are galaxies that cannot give good fits us- 
ing Newtonian dynamics plus dark matter (NGC 2976 and NGC 
7793) either, see de Blok et al. (2008). For the remaining galaxy 
(NGC 3198) there is indeed some tension between observations 
and the MOND fit, that might be explained by the presence of 
non-circular motions, a small distance (see Pig. |7j, or a value of 
ao at the lower end of our best-fit interval (see Fig. IS). In any 
case, further observations (constraining the distance) and mod- 
elling of NGC 3198 in the MOND context should thus lead to 
a benchmark test for MOND as a modification of gravity. But 
we also show that MOND, as an empirical law encompassing 
the behavior of the gravitational field on galaxy scales, whatever 
its cause, is still very successful and summarizes old and new 
scaling relations with a remarkable consistency (Sect. 4.5). 

We also conclude that, both from arguments of best-fit stellar 
mass-to-light ratios (Sect 4.2) and best-fit distances (Sect. 4.3), 
the simple yu-function is preferred over the standard one. As 
noted by Famaey & Binney (2005) and McGaugh (2008), this 
is also the case when fitting the terminal velocity curve of our 
own Milky Way galaxy. Angus, Famaey & Diaferio (2010) also 
reached the same conclusion from using temperature profiles of 
the X-ray emitting gas of a sample of clusters, and from as- 
suming that dark matter in MONDian galaxy clusters is made 
of 1 leV fermionic particle^] 

Let us however note that, in order to constrain /j. from 
galaxy rotation curves within the modified gravity framework 
of MOND (see Sect. 2), one should actually calculate predic- 
tions of the modified Poisson formulations of Bekenstein & 
Milgrom (1984) or Milgrom (2010) numerically for each galaxy 
model, and for each choice of parameters. Our present conclu- 
sion for THINGS galaxy rotation curves does hold only for the 



3 Let us note that such a hot dark matter component could also play 
a role in the strong and weak gravitational lensing of elliptical galaxies 



baryonic gravitational acceleration, the scatter is much lower (Ferreras et al. 2009, Tian et al. 2009) 
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Fig. 9. Mass discrepancy-acceleration relation using the rotation 
curve data of our sample. The gravitational acceleration gener- 
ated by baryons (g N ) is measured in km 2 s~ 2 kpc~' and they 
result from the fits made with the simple interpolation function 
and the distance constrained in Table 2. Black (open) circles rep- 
resent the data points with an uncertainty larger than 5%. The 
data points with an uncertainty smaller than 5% are shown as 
red (full) circles. 
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Fig. 10. Gravitational acceleration generated by phantom dark 
matter (gphantomDM, measured in km 2 s~ 2 pc~') versus radius (in 
pc). The values of gphantomDM result from the fits made with 
the simple interpolation function and the distance constrained 
in Table 2. 

modified inertia formulation for circular orbits given here by 
Eq. 1. 
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Table 2. Mass modelling results. M/Lj isk and M/L bu i 5C are in the 3.6/um band. The values called / are the ratio between the best-fit value and the 
"nominal" one (stellar M/L from the colours and the independently determined distance). The distance is in Mpc and^ 2 ^ is the reduced^ 2 (this 
number is useful to compare different models of the same galaxy but should not be used as a true probability indicator, see footnote[TJ. The values 
refer to the fits made with ao = 1.22 x 10~ 8 cm s~ 2 . 
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0.00 


8.90 




0.63 




// simnlp d mnstr 


1.71 


2.80 


0.00 


0.00 


9.56 


1.07 


0.58 




// simnlp ^ frpp 

LA O 1 J 1 J I J 1 v. (. J 1 1 H 


1.71 


2.80 


0.00 


0.00 


9.56 


1.07 


0.58 


»\ vJv. JUJJ 


»/ ctonn /V fi v a/1 
yU SLallU. u 11ACU 




^4 


/t6 


o 

U. JJ 


10 10 
1U. 1U 




? 6^ 
Z.O J 




jU SLallU. Ct COI1SLF. 


\J. 1 J 


CK 

U.7J 


o 

U.J J 


Al 
U.t-Z 


7 07 
/ .U / 


70 

u. / u 


Q7 




LL o LdllLI. U 11 ll 


84 


1 06 

1 .l/VJ 


57 


44 


6 55 


65 


91 
yj.y i 




/i simple fixed 


0.30 


0.38 


0.43 


0.33 


10.10 




1.80 




// simnlp /V ronstr 

Ll ..i 1 J 1 1 I J 1 v. (. J IW 1 1 1 1 L± > 


0.55 


0.70 


0.56 


0.43 


7.07 


0.70 


0.86 




/i simple free 


0.66 


0.84 


0.61 


0.47 


6.27 


0.62 


0.80 


\r rip ^ 1 QR 

IN. Vj^ J 1 70 


fj, SLallU. Ll HAeU 


n /to 


61 
W.u 1 






1 ^ 80 
1 J.oU 




6 97 
o.z / 




// ctonrl /V rr*nctr 
LL dLallU. Ll CUllftLl. 


67 


Rzt 






1? ^0 
1Z. JU 


8Q 
VJ.oy 


^ 81 
J.o 1 




// QtanH d i~v&& 

LL o LailLI. U 11 tl 


1 43 

1 ."-TO 


1 79 

i . / y 






8 38 

O.JO 


61 

l/.VJ 1 


1 56 




// simnlp d fixprl 

LA, O 1 J 1 J I J 1 1 (. J lli\l U 


0.37 


0.46 






13.80 




6.18 




// simnlp d ronstr 

LA. ^>1J11I./11 XA, V— V ' 1 I.I L± . 


0.48 


0.60 






12.30 


0.89 


3.79 




// simnlp d frpp 

LA ^>lJllL'il L,t 1 1 I I 


1.01 


1.26 






8.60 


0.62 


1.32 


Mnp ^6?1 


// ctofln /V flVPn 
LL SLallU. Ll 11ACU 


W.J 1 


86 
u.ou 






6 6/1 




70 
u. / u 




fj, SLallU. Ll COIlSLr. 


n fin 


1 0? 
1 .wz 






6 1 A 

O. It- 


Q9 

u.vz 


so 

U.JU 




// ct^inri d TYf*f* 

LL MtlllLl. Ll 11CC 


U.UU 


1 0? 






6 14. 


Q? 


SO 

U.JU 




yU simple d fixed 


0.37 


0.63 






6.64 




0.78 




/i simple d constr. 


0.44 


0.75 






6.11 


0.92 


0.55 




yU simple d free 


0.44 


0.75 






6.11 


0.92 


0.55 


NGC 2403 


yU stand, d fixed 


0.74 


1.80 






3.47 




2.12 




yU stand, d constr. 


0.62 


1.51 






3.76 


1.08 


1.43 




yU stand, ii free 


0.35 


0.85 






4.71 


1.46 


0.54 




yU simple ii fixed 


0.53 


1.29 






3.47 




2.29 




/i simple ii constr. 


0.45 


1.10 






3.76 


1.08 


1.73 




yU simple d free 


0.26 


0.63 






4.69 


1.46 


0.56 


NGC 7793 


yii stand, d fixed 


0.46 


1.48 






3.91 




6.21 



4 In NGC 2903 the best-fit M/L of the bulge is zero, but values larger than the best-fit M/L of the disk give almost equally good fits. 
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Table 2. continued. 



Name 


fit 
lit 


A/ill 


/disk ^/-Mmlge /bulge 


dist 


Jdist 


V 2 
'I red 




/i stand, d constr. 


0.39 


1.26 


4.30 


1.10 


5.40 




/i stand, d free 


0.14 


0.45 


7.02 


1.80 


3.08 




/i simple c/ fixed 


0.33 


1 .06 


3.91 




5.94 




/i simple <i constr. 


0.28 


0.90 


4.30 


1.10 


5.17 




/i simple c/ free 


U.1Z 


U.j9 


6.56 


1.6o 


j.UU 


NGC 2976 


/i stand, d fixed 


0.33 


0.60 


3.56 




1.73 




/i stand, d constr. 


0.28 


0.51 


3.90 


1.09 


1.65 




4. J J J? 

/i stand, a free 


0.04 


0.07 


9.30 


2.61 


0.70 




simple a fixed 


0.23 


0.42 


3.56 




1 £ 1 

1.61 




/i simple <i constr. 


0.20 


0.36 


3.92 


1.10 


1 CA 

1.50 




/i simple d free 




n no 


7 79 

/ . /Z 


9 17 


n 77 
u. / / 


DDO 154 


/i stand, d fixed 


0.00 


0.00 


4.30 




4.60 




yu stand. <f constr. 


0.77 


2.41 


3.23 


0.75 


0.36 




fi stand. <f free 


0.85 


2.66 


3.18 


0.74 


0.35 




/i simple <f fixed 


0.00 


0.00 


4.30 




6.73 




/i simple d constr. 


0.50 


1.56 


3.23 


0.75 


0.41 




yii simple d free 


0.79 


2.47 


3.04 


0.71 


0.33 
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